# library
library(tidyverse)
library(ggplot2)
library(dplyr)
library(ggbeeswarm)Tulai
# raed the data
tl <- read_csv("Tulai_Lithic_Assamblage.csv")# create a new data frame for plotting
tl1 <- tl %>%
select(AREA,
`Width (mm)`,
`Thickness (mm)`,
`Length (mm)`) tl1 %>%
filter(AREA == "A") %>%
select(AREA, `Width (mm)`) %>%
ggplot() +
aes(x = `Width (mm)`) +
geom_histogram(binwidth = 2) +
labs(x = "Width of Debitage (mm)",
y = "Frequency",
title = "Histogram of Debitage Width for Area A" )+
theme_minimal()#| label: Width of Artifacts (Area A)
#| fig-cap: "Width of Artifacts (Area A)"
#| fig-height: 5
tl1 %>%
filter(AREA == "TP1") %>%
select(AREA, `Width (mm)`) %>%
ggplot() +
aes(x = `Width (mm)`) +
geom_histogram(binwidth = 2) +
labs(x = "Width of Debitage (mm)",
y = "Frequency",
title = "Histogram of Width for Area TP1" )+
theme_minimal()tl1 %>%
filter(AREA == "D1, D3, D4, D5, C, B1, DH") %>%
select(AREA, `Width (mm)`) %>%
ggplot() +
aes(x = `Width (mm)`) +
geom_histogram(binwidth = 5) +
labs(x = "Width of Debitage (mm)",
y = "Frequency",
title = "Histogram of Width for other Areas" )+
theme_minimal()ggplot(data = tl1) +
aes(x = `Width (mm)`,
y = AREA) +
geom_boxplot() +
labs(x = "Width",
y = "Areas") +
theme_classic(base_size = 14)